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Abstract 

A measurement of the production cross-section for top quark pairs (it) in pp collisions at -\fs — 7 TeV is presented using data 
recorded with the ATLAS detector at the Large Hadron Collider. Events are selected in the single lepton topology by requiring 
\«J "an electron or muon, large missing transverse momentum and at least three jets. With a data sample of 35 pb -1 , two different 
2T» multivariate methods, one of which uses Z?-quark jet identification while the other does not, use kinematic variables to obtain cross- 
section measurements of erg — 187 ± ll(stat.)_j°(syst.) ± 6(lumi.) pb and cr„- = 173 ± 17(stat.)_Jg(syst.) ± 6(lumi.) pb respectively. 
O .The two measurements are in agreement with each other and with QCD calculations. The first measurement has a better a priori 
sensitivity and constitutes the main result of this Letter. 
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Measurements of the production and decay properties of top 
quarks are of central importance to the Large Hadron Collider 
(LHC) physics programme. Uncertainties on the theoretical 
predictions for the top quark pair production cross-section are 
now less than 10%, and comparisons with experimental mea- 
surements allow a precision test of the predictions of Quantum 
Chromodynamics. Furthermore, top quark pair production is an 
important background in many searches for physics beyond the 
Standard Model (SM). New physics may also give rise to ad- 
ditional it production mechanisms or modifications of the top 
quark decay channels, which can affect the measured tt cross- 
section. 

In the SM the tt production cross-section in pp collisions is 
calculated to be 165 ^[g pb JH S [H at a centre-of-mass en- 
ergy yfs — 7 TeV, assuming a top quark mass of 172.5 GeV. 
Top quarks are predicted to decay to a W-boson and a fo-quark 
(t — > Wb) nearly 100% of the time. Events with a tt pair can 
be classified as 'single lepton', 'dilepton', or 'all hadronic' ac- 
cording to the decays of the two VK-bosons: each can decay 
into quark-antiquark pairs (W — > qiq~2) or a lepton-neutrino pair 
(W — * £v). Events in the single lepton channel, when the lep- 
ton is an electron or a muon, are characterised by an isolated, 
prompt, energetic lepton, jets, and missing transverse momen- 
tum from the neutrino. At the Tevatron the tt cross-sections 
at sis = 1.8 TeV and at V* = 1 -96 TeV have been measured 
by CDF Jl H and D0 Jl 0] in most channels. ATLAS and 
CMS have measured the tt cross-section at -\fs — 7 TeV at the 
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This Letter describes measurements of the it cross-section in 
the single lepton plus jets channel with 35 pb 1 of data recorded 
by ATLAS in 2010. Taking advantage of the increased data 
sample, the measurement techniques developed in Ref. [ 8] were 
extended to employ kinematic likelihood discriminants to sep- 
arate signal from background and measure the cross-section. 
Two multivariate methods, one that includes fr-quark jet iden- 
tification (^-tagging) and one which does not, use several vari- 
ables each to discriminate it events from the background. The 
two analyses are sensitive to different sources of systematic un- 
certainty. For instance, the analysis without ^-tagging is more 
sensitive to the multijet background, whereas the analysis with 
b-tagging is sensitive to the background from W-boson produc- 
tion in association with b- and c-quarks. The clearer separation 
of signal and background leads to a smaller statistical uncer- 
tainty for the analysis with ^-tagging. Another significant dif- 
ference between the two measurements is that the analysis with 
^-tagging uses a profile likelihood that implements an in situ 
fit of the dominant systematic uncertainties, which improves its 
performance considerably. 



2. The ATLAS detector 

The ATLAS detector ll 211 consists of an inner tracking sys- 
tem (inner detector, or ID) surrounded by a thin superconduct- 
ing solenoid providing a 2 T magnetic field, electromagnetic 
and hadronic calorimeters and a muon spectrometer (MS). The 
ID consists of silicon pixel and microstrip detectors, surrounded 
by a transition radiation tracker. The electromagnetic calorime- 
ter is a lead/liquid-argon (LAr) detector. Hadron calorimetry is 
based on two different detector technologies, with scintillator 
tiles or LAr as active media, and with either steel, copper, or 
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tungsten as the absorber material. The MS includes three large 
superconducting toroids arranged with an eight-fold azimuthal 
coil symmetry around the calorimeters, and a system of three 
stations of chambers for the trigger and for track measurements. 

A three-level trigger system is used to select interesting 
events. The level- 1 trigger is implemented in hardware and 
uses a subset of detector information to reduce the event rate 
to a design value of at most 75 kHz. This is followed by two 
software-based trigger levels, level-2 and the event filter, which 
together reduce the event rate to about 200 Hz which is recorded 
for analysis. 

The nominal pp interaction point at the centre of the detec- 
tor is defined as the origin of a right-handed coordinate system. 
The positive x-axis is defined by the direction from the interac- 
tion point to the centre of the LHC ring, with the positive y-axis 
pointing upwards, while the z-axis is along the beam direction. 
The azimuthal angle (p is measured around the beam axis and 
the polar angle 9 is the angle from the z-axis. The pseudorapid- 
ity is defined as r\ — - In tan(#/2). 

3. Simulated event samples 

Monte Carlo (MC) simulation was used for various aspects of 
the analysis. The simulation consists of an event generator in- 
terfaced to a parton shower and hadronisation model, the results 
of which are passed through a full simulation of the ATLAS de- 
tector and trigger system lfl3l [3l ■ MC simulation was used 
when data-driven techniques were not available or to evaluate 
relatively small backgrounds and certain sources of systematic 
uncertainty. 

For the calculation of the acceptance of the it signal the next- 
to-leading order (NLO) generator MC@NLO v3.41 [ 15] was 
used with the top quark mass set to 172.5 GeV and with the 
NLO parton density function (PDF) set CTEQ66 0]. 

W- and Z-boson production in association with jets was sim- 
ulated with Alpgen v2.13, which implements the exact leading 
order (LO) matrix elements for final states with up to six partons 
and uses the 'MLM' matching procedure to remove the overlaps 
between samples with n and n + 1 final state partons [17]. The 
LO PDF set CTEQ6L1 [16] was used to generate W+jets and 
Z+jets events with up to five partons. Diboson, WW, WZ and 
ZZ events were generated with Herwig Il8lll9ll . Like the dibo- 
son production, single-top is also a relatively small background 
and is simulated using MC@NLO, invoking the 'diagram re- 
moval scheme' J20I1 to remove overlaps between single-top and 
tt final states. 

Unless otherwise noted, all events were hadronised with 
Herwig, using Jimmy 112 ill for the underlying event model. De- 
tails of the generator and underlying event tunes used are given 
in Ref. iQ. 



3.1. Systematic uncertainties on signal and background mod- 
elling 

The use of simulated tt samples to calculate the signal accep- 
tance gives rise to various sources of systematic uncertainty. 
These arise from the choice of the event generator and PDF 



set, and from the modelling of initial and final state radiation 
(ISR and FSR). The uncertainties due to the choice of genera- 
tor and parton shower model were evaluated by comparing the 
results obtained with MC@NLO to those of Powheg H23I1 . with 
events hadronised with either Herwig or Pythia 1240. The un- 
certainty due to the modelling of ISR/FSR was evaluated using 
the AcerMC generator B25I1 interfaced to Pythia and by varying 
the parameters controlling the ISR/FSR emission by a factor of 
two up and down. The variation ranges used are comparable 
to those in M26I1 for ISR and 12711 for FSR. Finally, the uncer- 
tainty in the PDF set used to generate tt samples was evaluated 
using a ran ge o f current PDF sets with the procedure described 
inRef. HIM 113] • 

The production of the W+jets background based on MC sim- 
ulation has uncertainties on the total cross-section, on the con- 
tribution of events with jets from heavy-flavour (b, c) quarks, 
and on the shape of kinematic distributions. The predictions 
of the total cross-section have uncertainties of order 50% [31], 
increasing with jet multiplicity. Total W+jets cross-section pre- 
dictions were not used in the cross-section measurement as 
this background was extracted from the fit to the data (see 
Section [7]), but were used in the MC simulation shown in 
Figs. Q] to |U A combination of the fitting method described in 
[32] and a counting method described here, both relying upon 
final states with one and two jets, was used to estimate the heavy 
flavour fractions in W+jets events. Since these bins are domi- 
nated by W+jet events, the total W+jet contribution to these 
events can be obtained, both with and without requiring at least 
one ^-tagged jet. These four numbers are then used to constrain 
the following four event types which make up the W+jets sam- 
ple: W + bb, W+cc, W+c and W+light flavours. Additionally it 
was assumed that the k-factors for W + bb and W + cc are equal. 
MC simulation with Alpgen was used to estimate the ^-tagging 
efficiencies for each sub-sample as well as to extrapolate from 
the one-jet to the two-jet bin. The dominant uncertainties in this 
method arise from jet energy scale and ^-tagging uncertainties. 
As a result of this study, it was found that the W + bb and W + cc 
sub-samples of events in the Alpgen MC simulation were to be 
rescaled by 1.30±0.65, whereas W + c events were rescaled by 
1 .0±0.4. An additional 25% relative uncertainty per jet bin was 
assigned to these flavour fractions when applied to the signal 
region based upon studies with Alpgen MC simulation. 

The uncertainty on the shape of W+jets kinematic distribu- 
tions was assessed by changing the factorisation and renormal- 
isation scales by a factor of two up and down; and by varying 
the minimum pj of the final state quarks and gluons from 10 to 
25 GeV, with 15 GeV being the default. 

For the smaller backgrounds arising from Z+jets, single-top 
and diboson production, only the overall normalisation uncer- 
tainties were considered, taken to be 30% for Z+jets produc- 
tion, 10% for single-top production, determined from compar- 
isons of MCFM Q and MC@NLO predictions, and 5% for 
diboson production, determined from MCFM studies of scale 
and PDF uncertainties. 
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4. Object selection 

Single lepton ti events are characterised by the presence 
of an electron or muon, jets, and missing transverse momen- 
tum, which is an indicator of undetected neutrinos, in the fi- 
nal state. The events used in this analysis were triggered by 
single-lepton triggers. The electron trigger required a level- 1 
electromagnetic cluster in the calorimeter with transverse mo- 
mentum Ej > 10 GeV. A more refined cluster selection was 
applied in the level-2 trigger, and a match between the electro- 
magnetic cluster and an ID track was required in the event filter. 
The muon trigger required a track with transverse momentum 
pr > 10 GeV in the muon trigger chambers at level- 1, matched 
to a muon of pj > 13 GeV reconstructed in the precision cham- 
bers and combined with an ID track at the event filter. 

The same object definition used for the previous it cross- 
section measurement |@1 was used in this analysis, except for 
more stringent electron selection criteria and ID track quality 
requirements for muons. Electron candidates were defined as 
electromagnetic clusters consistent with the energy deposition 
of an electron in the calorimeters and with an associated well- 
measured track. They were required to satisfy pj > 20 GeV 
and |?7ciusteil < 2.47, where fluster is the pseudorapidity of the 
cluster associated with the candidate. Candidates in the barrel 
to endcap calorimeter transition region 1.37 < |77ciusterl < 1-52 
were excluded. Muon candidate tracks were reconstructed from 
track segments in the different layers of the muon chambers. 
These segments were combined starting from the outermost 
layer, with a procedure that takes material effects into account, 
and matched with tracks found in the inner detector. The final 
candidates were refitted using the complete track information 
from both detector systems and required to satisfy pj > 20 
GeV and \rj\ < 2.5. 

To further reduce background from leptons produced in 
heavy-flavour or in-flight hadron decays the selected leptons 
were required to be 'isolated'. For electrons the transverse mo- 
mentum, Ej, deposited in the calorimeter cells inside an iso- 
lation cone of size AR = V(A?7) 2 + (A0) 2 = 0.2 around the 
electron position was corrected to take into account the leakage 
of the electron energy into this cone. The remaining Ej was 
required to be less than 4 GeV. Muons were required to have a 
distance AR greater than 0.4 from any jet with pj > 20 GeV, 
which suppresses muons from heavy-flavour decays inside jets. 
Furthermore, the calorimeter transverse momentum in a cone 
of size AR = 0.3 around the muon direction was required to be 
less than 4 GeV, and the sum of track transverse momenta, other 
than the muon track, in a cone of size AR = 0.3 was required to 
be less than 4 GeV. 

Pure samples of prompt muons and electrons were obtained 
from Z-boson events in the data and were used to correct the 
lepton trigger, and the reconstruction and selection efficiencies 
in MC simulation to match those in the data. The corrections 
were found to be small. 

Jets were reconstructed ll34ll with the anti-fc, algorithm [35, 
l36ll with radius parameter 0.4 from clusters of adjacent 
calorimeter cells. If the closest object to an electron candi- 
date (before the above electron isolation requirement) was a jet 



within a distance AR < 0.2, the jet was removed. The jet en- 
ergy scale (JES) and its uncertainty were derived by combining 
information from test-beam data, LHC collision data and simu- 
lation. The JES uncertainty was found to vary from 2% to 7% 
as a function of jet p T and 77 113711 . 

Jets arising from the hadronisation of Z?-quarks were identi- 
fied using an algorithm (JetProb) l38ll which relies upon the 
transverse impact parameter do of each track in the jet: this is 
the distance of closest approach in the transverse x-y plane of 
a track to the primary vertex. It is signed with respect to the 
jet direction: the sign is positive if the track crosses the jet axis 
in front of the primary vertex, negative otherwise. The signed 
impact parameter significance, do/(Td , of each selected track is 
compared to a resolution function for prompt tracks, to assess 
the probability that the track originates from the primary vertex. 
Here, <x rfo is the uncertainty on do. The individual track proba- 
bilities are then combined into a global probability that the jet 
originates from the primary vertex. The simulated data were 
smeared to reproduce the resolution found in collision data. 

The /^-tagging efficiencies and mistag rates were calibrated 
with data for a wide range of ^-tagging efficiency requirements. 
The efficiency was measured in a sample of jets containing 
muons, making use of the transverse momentum of the muon 
relative to the jet axis. The mistag rates were measured on an 
inclusive jet sample with two methods, one using the invariant 
mass spectrum of tracks associated to reconstructed secondary 
vertices to separate light- and heavy-flavour jets, and the other 
based on the fraction of secondary vertices in data with neg- 
ative decay-length significance. The results of these measure- 
ments were applied in the form of /^-dependent scale factors 
to correct the ^-tagging performance in simulation to match the 
data. For a ^-tagging efficiency around 50%, the scale factor 
was found to be approximately 0.9 in all bins of jet pj, and 
the relative ^-tagging efficiency uncertainty was found to range 
from 5% to 14% depending on the jet pj 03811 . The mistag rate 
and mistag scale factors are approximately 1% and 1.1, respec- 
tively, in the jet pj region of interest, 20 < pj < 100 GeV. The 
analysis including ^-tagging used the probabilities returned by 
the JetProb algorithm as a discriminating variable, as explained 
in Section [7] 

The reconstruction of the missing transverse momentum 
£™ ss [39] was based upon the vector sum of the transverse mo- 
menta of the reconstructed objects (electrons, muons, jets) as 
well as the transverse energy deposited in calorimeter cells not 
associated with these objects. The electrons, muons and jets 
were used in the £™ lss calculation consistently with the defini- 
tions and uncertainties stated above. 



5. Event selection 

Events that passed the trigger selection were required to 
contain exactly one reconstructed lepton with pj > 20 GeV, 
matching the corresponding event filter object. Selected events 
were required to have at least one reconstructed primary ver- 
tex with at least five tracks. Events were discarded if any jet 
with pj > 20 GeV was identified to be due to calorimeter noise 
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or activity out of time with respect to the LHC beam cross- 
ings. The £™ ss was required to be greater than 35 (20) GeV 
in the electron (muon) channel and the transverse mass con- 
structed from the lepton and £™ ss transverse momentum vec- 
tors was required to be greater than 25 GeV (60 GeV -Ej lss ) in 
the electron (muon) channel. The muon requirement is referred 
to as the 'triangular cut'. The requirements were stronger in 
the electron channel to suppress the larger multijet background. 
Finally, events were required to have three or more jets with 
Pt > 25 GeV and \rj\ < 2.5. The selected events were then clas- 
sified by the number of jets fulfilling these requirements and by 
the lepton flavour. Table[T]shows the number of selected events 
in the data in the electron and muon channels, together with the 
SM expectations for the signal and the different backgrounds. 
All predictions were obtained from MC simulation except the 
multijet background estimate which was obtained from data as 
described in the next section. 



Table 1: Number of observed events in the data in the electron and 
muon channels after the selection cuts as a function of the jet multiplic- 
ity. The expected signal and background contributions are also given. 
All simulated processes are normalized to theoretical SM predictions, 
except the multijet background which uses the normalisation presented 
in Sec. [6] The quoted uncertainties include statistical, systematic and 
theoretical components, except for the multijet background. All num- 
bers correspond to an integrated luminosity of 35 pb _1 . 



Electron channel 


3 jets 


4 jets 


> 5 jets 


tt 


117 ± 


16 


109 ± 


15 


76 ± 


19 


W+jets 


524 + 


225 


124 ± 


77 


35 + 


23 


Multijet 


64 ± 


32 


12 ± 


6 


8± 


4 


Single top 


21 ± 


5 


7 ± 


3 


3± 


2 


Z+jets 


60 ± 


28 


21 ± 


15 


8± 


6 


Diboson 


9± 


3 


1.9 ± 


1.5 


0.4 ± 


0.8 


Predicted 


795 ± 


236 


275 ± 


84 


130 + 


35 


Observed 


755 


261 


123 


Muon channel 


3 jets 


4 jets 


> 5 jets 


tt 


165 ± 


22 


156 ± 


18 


108 ± 


27 


W+jets 


976 ± 


414 


222 ± 


139 


58 + 


38 


Multijet 


79 ± 


24 


18 + 


6 


11 ± 


3 


Single top 


31 ± 


7 


10 ± 


4 


4± 


2 


Z+jets 


58 ± 


26 


14 ± 


10 


5± 


4 


Diboson 


16 ± 


4 


3± 


2 


0.6 ± 


0.8 


Predicted 


1325 ± 


422 


423 ± 


143 


186 ± 


51 


Observed 


1289 


436 


190 



6. Background evaluation 

The main backgrounds to tt signal events in the single lep- 
ton plus jets channel arise from W-boson production in associ- 
ation with jets, in which the W decays leptonically, and from 
multijet production. Smaller backgrounds arise from Z+jets, 
diboson and single-top production. These smaller backgrounds 
have been estimated from MC simulation and normalised to the 
latest theoretical predictions, as discussed in Section[3] 



The W+jets background is difficult to predict from theory, 
particularly in the high jet-multiplicity bins. A data-driven 
cross-check following methods similar to those described in 
Ref. [8] was therefore performed. The results obtained with 
data were found to agree with the MC predictions within the 
uncertainties. Both analyses presented here rely on the assump- 
tion that the MC simulation correctly describes the kinematic 
properties of the W+jets events, whereas the normalisation of 
the W+jets cross-section was fitted from the data, as described 
in Section|7] In the analysis using ^-tagging the theoretical un- 
certainty on the normalisation was used as a constraint in the 
fit, whereas in the other analysis it was allowed to vary freely. 

The multijet background was measured with a data-driven 
approach. In the muon channel, the background from multi- 
jet events is dominated by 'non-prompt' muons arising from 
the decay of heavy-flavour hadrons, in contrast to the tt signal 
where muons arise from the 'prompt' decays of W-bosons. The 
multijet background can be estimated by defining two samples 
of muons, 'loose' and 'tight'. The tight sample is the one de- 
fined in the event selection described above, whilst the loose 
sample satisfy the same criteria except the muon isolation re- 
quirements. Since the reconstructed muons from background 
are associated with jets, they tend to be much less isolated than 
the leptons in tt decays. Any sample of muons is composed of 
prompt and non-prompt muons and it is assumed that the tight 
muon sample is a subsample of the loose sample: 



xrloose _ xrloose , xrloose 

— prompt non-prompt' 



W tight = e t N loose + e t W loose 

Jy e prompt J » prompt ' e non-prompt J » non-prompt' 



(1) 



where A^°°!f prompt is the number of loose, non-prompt muons 
(with the other Ny's defined similarly) and e prom pt (^non-prompt) 
represents the probability for a prompt (non-prompt) muon that 
satisfies the loose criteria to also satisfy the tight ones. The 
probability e P rompt was measured from the data using high- 
purity samples dominated by Z-bosons decaying into muons. 
The probability e non -prompt for a non-isolated lepton to pass the 
isolation cuts was measured by defining control samples dom- 
inated by multijet events. Two different control samples were 
defined to have at least one jet plus a muon (i) with high impact 
parameter significance or (ii) with low transverse mass of the 
muon-£' mlss system plus reversed triangular cut. These control 
samples gave consistent results. Contamination of the multijet 
control samples by muons from W and Z events was determined 
from MC simulation. The results of these studies are e n on-prompt 
and e P rompt as a function of the muon 77, from which the multijet 
background expectations can be obtained as a function of any 
variable. A 30% systematic uncertainty was assigned to this 
estimate based on the observation that the method gives agree- 
ment to within 30% across the different jet multiplicities. 

In the electron channel, the multijet background also in- 
cludes photons inside jets undergoing conversions into electron- 
positron pairs and jets with high electromagnetic fractions. A 
different method was used, based on a binned likelihood fit of 
the £ miss distribution in the region £ miss < 35 GeV The data 
was fitted to the sum of four templates: multijet, tt, W+jets and 
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Z+jets. The templates for the latter three processes were ob- 
tained from MC simulation whereas the multijet template was 
obtained from the data in a control region defined by the full 
event selection criteria except that the electron candidate fails 
one or more of the identification cuts. The multijet background 
was obtained by extrapolating the fraction of multijet events 
from the fit at low £™ ss to the signal region at high £™ ss . Sev- 
eral choices of electron identification cuts were considered and 
the largest relative uncertainty among these (50%) was used as a 
conservative estimate of the systematic uncertainty of this back- 
ground evaluation. 

7. Cross-section extraction 

The tt production cross-section was extracted by exploiting 
the kinematical properties of tt events compared to those from 
the dominant background (W+jets) by means of likelihood dis- 
criminants (D) constructed from several variables. Templates 
of the distributions D for signal and all background samples 
were created using the TMVA package [40]. The variables were 
selected for their good discriminating power, small correlation 
with each other, and low sensitivity to potentially large uncer- 
tainties such as jet energy calibration. The variables are: 

• The pseudorapidity rj of the lepton, since leptons produced 
in tt events are more central than those in W+jet events. 

• The aplanarity M, defined as 3/2 times the small- 
est eigenvalue of the momentum tensor M i; = 
£ d^tt, p^p^lY^jf"* p 2 k , where pik is the z'-th momentum 
component of the k-th object and pi is the modulus of its 
momentum. The lepton and the four leading jets are the 
objects included in the sum. To increase the separation 
power of the aplanarity distribution, the transformed vari- 
able exp(-8 x was used. This variable exploits the 
fact that tt events are more isotropic than W+jets events. 

• The charge of the lepton gi ep ton, which uses the fact that a 
sample of tt events should contain the same number of pos- 
itively and negatively charged leptons, while W+jet events 
produce an excess of positively charged leptons in pp col- 
lisions. 

• Hjt, p , defined as the sum of the transverse energies of 
the third and fourth leading jets normalised to the sum 
of the absolute values of the longitudinal momenta of the 
four leading jets, the lepton and the neutrino, Hjj p = 

Z£=3 \Pti\/ HijST' \Pz,j\> where pj is the transverse mo- 
mentum and p z the longitudinal momentum. The longi- 
tudinal momentum of the neutrino was obtained using the 
quadratic W mass constraint and taking the solution with 
the smaller neutrino p z value. To increase the separation 
power of the Hjj p distribution, the transformed variable 
exp (-4 x Hjj p ) was used. 

• The average wjp of wjp = - log 10 P/ for the two jets with 
lowest Pi in the event. P/ is the probability for a jet to 
be a light jet from the JetProb ^-tagging algorithm. These 



correspond to the jets that have the highest probability to 
be heavy-flavour jets. 

Two complementary analyses were performed, one which re- 
lied upon the use of Z?-tagging information (i.e. the variable vvjp) 
and one which did not. We refer to the analyses as 'tagged' and 
'untagged', respectively. The untagged analysis employed the 
first three variables, whereas the tagged analysis did not con- 
sider the lepton charge but used Hj t, p and vvjp. wjp was not 
included in the three-jet bin. Figures Q] to [4] show the distribu- 
tions of the discriminating variables for the selected data super- 
imposed on the signal and background SM predictions for the 
different jet multiplicities. 

The tt cross-section was extracted by means of a likelihood 
fit of the signal and background discriminant distributions to 
those of the data. The fit yields the fractions of tt signal and 
backgrounds in the data sample. The fit was performed simul- 
taneously to four samples (three-jet exclusive and four-jet inclu- 
sive, electron and muon) in the untagged analysis and six sam- 
ples (three-jet exclusive, four-jet exclusive and five-jet inclu- 
sive, electron and muon) in the tagged analysis, as these were 
the combinations that provided maximum sensitivity. The dis- 
criminants were built separately for each jet multiplicity and 
lepton flavour subsample, and the different channels were com- 
bined in the likelihood fit by multiplying the individual likeli- 
hood functions. 

The normalisation of the tt signal templates is the parame- 
ter of interest in the fit and was allowed to vary freely in both 
analyses. The tt cross-section was assumed to be common to 
all channels and the number of tt events in each subsample re- 
turned by the fit was related to the tt cross-section by the ex- 
pression aft = N S i g / (J Jjdt X e S igj, where N s i g is the number of 

tt events, J Hdt is the integrated luminosity and e slf is the prod- 
uct of the signal acceptance, selection efficiency and branching 
ratio, obtained from tt simulation. The normalisation of the 
backgrounds was treated differently in the two analyses. In the 
untagged analysis the multijet and small backgrounds (single- 
top, diboson and Z+jets production) were fixed in the fit to their 
expected contributions, whereas the W+jets background was al- 
lowed to vary freely in each channel. In the tagged analysis 
all backgrounds were allowed to vary within the uncertainties 
of their assumed cross-sections, described in Sections [3] and [6] 
These uncertainties were used as Gaussian constraints on the 
cross-section normalisation. The robustness of this fitting ap- 
proach was checked with ensemble tests. The central value and 
uncertainties returned by the fit were shown to be unbiased for 
a wide range of input cross-sections. 

8. Systematic uncertainties 

The evaluation of the systematic uncertainties was per- 
formed differently in the two analyses. The untagged analysis 
performed pseudo-experiments (PEs) with simulated samples 
which included the various sources of uncertainty. For exam- 
ple, for the JES uncertainty, PEs were performed with jet ener- 
gies scaled up and down according to their uncertainties and the 
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Figure 1 : Input variables to the likelihood discriminants in the exclusive three- 
jet bin for the muon channel: lepton r\ (top), exp(-8 X l7{) (second from top), 
lepton charge (third from top) and exp(-4x//x,3p) (bottom). All simulated pro- 
cesses are normalized to theoretical SM predictions, except the multijet back- 
ground which uses the normalisation presented in Sec. [6] The two top distribu- 
tions are used in the untagged and the tagged analyses, the third distribution in 
the untagged analysis, and the bottom distribution in the tagged analysis. 



Figure 2: Input variables to the likelihood discriminants in the inclusive four-jet 
bin for the electron channel: lepton r\ (top), exp(-8 X LK) (middle) and lepton 
charge (bottom). All simulated processes are normalized to theoretical SM 
predictions, except the multijet background which uses the normalisation pre- 
sented in Sec. [6] These distributions are used in the untagged analysis. 
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Figure 3: Input variables to the likelihood discriminants in the exclusive four- Figure 4: Input variables to the likelihood discriminants in the inclusive five-jet 

jet bin for the muon channel: lepton r\ (top), exp(-8 X J?() (second from top), bin for the electron channel: lepton r\ (top), exp(-8 X J?() (second from top), 

exp(-4 X H~T3 p ) (third from top) and wjp (bottom). All simulated processes exp(-4 X //t,3/j) (third from top) and wjp (bottom). All simulated processes 

are normalized to theoretical SM predictions, except the multijet background are normalized to theoretical SM predictions, except the multijet background 

which uses the normalisation presented in Sec. [6] These distributions are used j which uses the normalisation presented in Sec. [6] These distributions are used 

in the tagged analysis. in the tagged analysis. 



impact on the cross-section was evaluated. The tagged analysis, 
on the other hand, accounted for most of the changes in the nor- 
malisation and shape of the templates due to systematic uncer- 
tainties by adding 'nuisance' terms to the fit (4JJ. Templates of 
the samples with one standard deviation 'up' and 'down' vari- 
ations of the systematic uncertainty source under study were 
generated in addition to the nominal templates. The fit interpo- 
lated between these templates with a continuous parameter by 
means of a Gaussian constraint. Before the fit, the constraint 
was such that the mean value was zero and the width was one; a 
fitted width less than one means that the data were able to con- 
strain that particular source of uncertainty. The effects due to 
the modelling of the W+jets and multijet background shapes, 
initial and final state radiation, parton density function of the 
tt signal, NLO generator, hadronisation and template statistics 
cannot be fully described by a simple linear parameter control- 
ling the template interpolation. As a consequence, they were 
not treated as nuisance terms but obtained by performing PEs 
with modified simulated samples, as was done in the untagged 
analysis. 

The nuisance parameters of the systematic uncertainties were 
all fitted together taking into account the correlations among 
them in the minimisation process. As a consequence, the un- 
certainties on the fitted quantities obtained from the fit include 
both the statistical and the total systematic components. There- 
fore, to obtain an estimation of the individual contributions to 
the total uncertainty in the tagged analysis, each individual sys- 
tematic uncertainty was obtained as the difference in quadrature 
between the total uncertainty and the uncertainty obtained after 
having fixed the corresponding nuisance parameter to its fitted 
value. The central values of the nuisance parameters after the fit 
agreed with their input values. The fit was cross-checked using 
PEs where the starting value of the nuisance parameters was 
different than the nominal value. The result was found to be 
unbiased. In addition, large variations of the kinematic depen- 
dence of the nuisance parameters {e.g. the JES as a function of 
the jet pj) were considered and resulted in a negligible impact 
on the result of the fit. 

The systematic uncertainties on the cross-section for both 
methods are summarised in Table|2] The dominant effects in the 
untagged analysis were JES, multijet and W+jets backgrounds 
shape and ISR/FSR. The latter was also important in the tagged 
analysis, together with the uncertainty related to the signal MC 
generator. In addition, this analysis was sensitive to effects re- 
lated to ^-tagging, specifically the determination of the heavy- 
flavour content of the W+jets background and the calibration of 
the ^-tagging algorithm itself. The luminosity uncertainty was 
3.4% [42, Si. 



Several cross-checks of the cross-section measurements were 
performed. These included the results of the likelihoods applied 
to individual lepton channels and tt cross-section measurements 
done with simpler and complementary approaches, including 
cut-and-count methods and fits to kinematic variables such as 
the reconstructed top mass. These cross-checks gave consistent 
results within the uncertainties. 



Table 2: Statistical and systematic uncertainties on the measured tt 
cross-section in the untagged and tagged analyses. Multijet and small 
backgrounds normalisation uncertainties are already included in the 
statistical uncertainty (a/i) in the tagged analysis. W+jets heavy- 
flavour content and i-tagging calibration do not apply (n/a) to the un- 
tagged analysis. The luminosity uncertainty is not included in the ta- 
ble. 



Method 


Untagj 


;ed 


Ta 


gged 


Statistical Error (%) 


+ 10.1 


-10.1 


+5.8 




5.7 


Object selection (%) 


JES and jet energy resolution 


+4.1 


-5.4 


+3.9 




2.9 


Lepton reconstruction, 












identification and trigger 


+ 1.7 


-1.6 


+2.1 




1.8 


Background modelling (%) 


Multijet shape 


+3.5 


-3.5 


+0.8 




•0.8 


Multijet normalisation 


+ 1.1 


-1.2 




a/i 




Small backgrounds norm. 


+0.6 


-0.6 




a/i 




W+jets shape 


+3.9 


-3.9 


+ 1.0 




1.0 


W+jets heavy-flavour content 


n/a 




+2.7 




•2.4 


^-tagging calibration 


n/a 




+4.1 




3.8 


signal modelling (%) 


ISR/FSR 


+6.3 


-2.1 


+5.2 




5.2 


NLO generator 


+3.3 


-3.3 


+4.2 




4.2 


Hadronisation 


+2.1 


-2.1 


+0.4 




0.4 


PDF 


+ 1.8 


-1.8 


+ 1.5 




1.5 


Others (%) 


Simulation of pile-up 


+ 1.2 


-1.2 


< 


0.1 




Template statistics 


+ 1.3 


-1.3 


+ 1.1 




1.1 


Systematic Error (%) 


+ 10.5 


-9.4 


+9.7 




9.0 



9. Results and conclusions 

The results of the likelihood fits applied to the data are 
shown in Figs. [5] and [6] where the distributions of the dis- 
criminants in data are overlaid on the fitted discriminant dis- 
tributions of the signal and backgrounds. The final measured 
cross-section results are: cr r j = 173 ± 17(stat.) ^(syst.) ± 
6(lumi.) pb = 173^4 pb in the untagged analysis and <x„- = 
187+ 1 l(stat.) +^(syst.)±6(lumi.) pb = 187^ pb in the tagged 
analysis. The two measurements are in agreement with each 
other. The latter has a better a priori sensitivity and thus con- 
stitutes the main result of this Letter. It is the most precise tt 
cross-section measurement at the LHC published to date and is 
in good agreement with the SM prediction calculated at NLO 
plus next-to-leading-log order 165 +}i pb lUlHI 
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Figure 5: Untagged analysis: (Top) The distribution of the likeli- 
hood discriminant for data superimposed on expectations for signal 
and backgrounds, scaled to the results of the fit. The left bins corre- 
spond to the muon channel and the right bins to the electron channel. 
(Bottom) The ratio of data to fit result. 



Figure 6: Tagged analysis: (Top) The distribution of the likelihood dis- 
criminant for data superimposed on expectations for signal and back- 
grounds, scaled to the results of the fit. The left bins correspond to the 
muon channel and the right bins to the electron channel. (Bottom) The 
ratio of data to fit result. 
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